FINE STRUCTURE OF THE ZEROS OF 
ORTHOGONAL POLYNOMIALS, 
IV. A PRIORI BOUNDS AND CLOCK BEHAVIOR 



YORAM AND BARRY SIMON^ ^ 

Dedicated to Percy Deift on the occasion of his sixtieth birthday 

Abstract. We prove locally uniform spacing for the zeros of or- 
thogonal polynomials on the real line under weak conditions (Ja- 
cobi parameters approach the free ones and are of bounded vari- 
ation). We prove that for ergodic discrete Schrodinger operators, 
Poisson behavior implies positive Lyapunov exponent. Both re- 
sults depend on a priori bounds on eigenvalue spacings for which 
we provide several proofs. 



1. Introduction 

Our primary goal in this paper concerns the fine structure of the 
zeros of orthogonal polynomials on the real line (OPRL), although we 
will say something about zeros of paraorthogonal polynomials on the 
unit circle (POPUC) (see Section[TUJ. Specifically, dfi will be a measure 
of compact support which is nontrivial (i.e., not supported on a finite 
set), usually a probability measure. Pn{x) or Pn{x]dfi) will be the 
monic orthogonal polynomials and Pn{x) = Pn/\\Pn\\ the orthonormal 
polynomials. The Jacobi parameters, {an,bn}'^=i, are defined by 

xPn{x) = Pn+l{x) + bn+lPn{x) + a'i^„-l(x) (1.1) 

= 0, 1, 2, . . . , where P^i{x) = 0. It follows that (when fi{R) = 1) 

||P„|| = ai . . .a„ (1.2) 

so 

xpn{x) = an+iPn+lix) + hn+iPn{x) + a„p„_i(a;) (1.3) 
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and X has a matrix representation in the orthonormal basis {pn}'^=o, 

/h «i ■■■\ 

J = Ui 62 a2 ■ ■ ■ (1.4) 

called the Jacobi matrix. The finite Jacobi matrix, Jn;F, is the n x n 
submatrix of J in the upper left corner. It is easy to see that (see (b) 
of the Appendix) 

Pn{x) = det(xl - Jn;F) (1-5) 

Let {x^J^^}^^^ be the zeros of Pn{x), so p.5|) says that the x^"^ are 
eigenvalues of Jn-p- Let dun be the pure point probability measure that 
gives weight 1/n to each x^""^ . We say the density of states exists if di^n 
has a weak limit duoo- By fll.5|l . one sees that 

/ x^di/„(x) = -Tr(4V) (1.6) 
J n 

which is often useful. 

The existence of the limit for a large class of regular measures on 
[—2,2] goes back to Erdos-Turan [9]j. Nevai [SI] realized all that was 
used is 

lim \bn\ + |a„ - 1| = (1.7) 



n— >oo 



(following the convention in the Schrodinger operator community, we 
use a„ = 1 as a "free" case, while the OP community uses a„ = |). 
Indeed, 

Theorem 1.1 (known). If d^i is a measure where p.7p holds, the den- 
sity of states exists and is given by 

du = 71-^ {VA-x^y\i^2,2] dx (1.8) 

The modern proof is not hard. For a„ = 1, 6„ = 0, the OPRL are 
explicitly given by 

smfc' 

from which one computes du exactly for this case. If Jo,n;F is the 
corresponding cutoff Jq, then ()1.7|) implies 

iTr(J„^^- jJ„^^)-.0 (1.9) 

for each k = 0, 1, 2, . . . which, by p.6|) . implies J x'^ dvn has the same 
limits as for the free case. 

Another case where it is known that dv exists is when 
samples of an ergodic family, that is, a^n^ = f{T'^u), b^^ = g{T^uj) 
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with T: Q Q an ergodic transformation on {Q,dp), a probability 
measure space. In that case, it is known (going back to the physics 
hterature and proven rigorously by Pastur [3H], Avron-Simon jl], and 
Kirsch-Martinelh |2H]): 

Theorem 1.2 (known). For ergodic Jacobi matrices, dvn,Lo has a limit 
dv for a.e. uo and dv is a.e. uj -independent. 

Again, the proof uses ()1.6|) plus in this case that, by ergodicity, 
^Tr(J^.^) has a limit a.e. by the Birkhoff ergodic theorem. 

The most important examples of the ergodic case are periodic, almost 
periodic, and random. 

One easily combines the two ideas to see that dv exists (and does 
not depend on 5a„, if 

an = al^'^ + 6an bn = b^^^ + 5hn 

with a^n\b^n'^ ergodic and |5a„| + |56„| — > 0. 

These results describe the bulk features of the zeros. Here we are 
interested in the fine structure, on the level of individual eigenvalues; 
specifically, the focus in jUl |12l US] and a main focus in this paper is 
what we call clock behavior, that the spacing locally is equal spacing. 
The term clock comes from the case of orthogonal polynomials on the 
unit circle (OPUC) where dv is typically Lebesgue measure on a circle 
and the equal space means the zeros look like the numbers on a clock. 

In order for the density of zeros to be dv, the equal spacing must 
be l/{dv/dE). The symmetric derivative dv/dE exists for a.e. E and, 
of course, (dv/dE) dE is the a.e. part of dv. It is known (see, e.g., 
Avron-Simon jT]) that dv has no pure points and, in many cases, it is 
known that dv/dE is a continuous function, or even C°°. To be formal, 
we define first 

Definition. Let Eq G supp (dv). We let zi^\Eo) be the zeros nearest 
Eq so that 

4-'H£^o) < zi-'\Eo) <E,< z^\E,) < z^:\E,) <... 

if such zeros exist. If Eq G [supp (rfz/ then z^^IEq) exists for j fixed 
and n large. 

Definition. Let Eq G supp{dv). We say there is weak clock behavior 
at Eq if dv/dEo exists and 

hm n[z^:\Eo) - zi'\Eo)] % {E,) = 1 (1.10) 
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We say there is strong clock behavior at Eq if du/dEo exists, p.lOj) 
holds, and for j = 1, ±2, ±3, . . . fixed, 

hm n[z^\E,) - z\i^'\E,)] ^ (Eq) = 1 (1.11) 

Definition. We say there is uniform clock behavior on [a, (3] if du/dE 
is continuous and nonvanishing on [a,j3] and 



lim 



sup 



E, E' are successive zeros of p„ in [a, P] 



(1.12) 

It is obvious that uniform clock behavior implies strong clock behav- 
ior at each interior point. 

In the earlier papers in this series that discussed clock behavior for 
OPRL im 031; there was a technical issue that severely limited the 
results in general situations. In all Jost function-type analysis 

was used to show that a suitably rescaled p„ converged, that is, 

CnPn (^Eo + n{x - Eo) ^ {Eo) + ^ fM (1.13) 

where < c/n and f^o has zeros at 0, ±1, ±2, .... This would naively 
seem to show that Pn has clock-spaced zeros and, indeed, it does imply 
at least one zero near Eo + ^n + f (^)~^ consistent with clock spacing. 
The snag involves uniqueness, for the function 

sin(7ra;) sin(n^7ra;) (1-14) 

n 

has a limit like foo but has more and more zeros near x = 0. That is. 



one needs to prove uniqueness of the zeros near Eq + C,n + |:(^)~^- 

In previous papers in this series, two methods were used to solve this 
uniqueness problem. One relied on some version of the argument prin- 
ciple, essentially Rouche's theorem. This requires analyticity which, 
typically, severely restricts what recursion coefficients are allowed. In 
the case of OPUC where one needs to control zeros in the complex 
plane, some kind of analyticity argument seems to be necessary. The 
second method relies on the fact that if (jl.l3|) also holds for derivatives 
and f'aoij) 7^ 0, then there is a unique zero. This argument also requires 
extra restrictions on the recursion coefficients, albeit not as severe as 
the analyticity requirement. For example, in |1T], one only needed 
Er=i|an - 1| + \bn\ < oo to get (ITT!^ for E e [-2 + e, 2 - e]. But to 
control derivatives, we needed J2'^=i ^(kn — 1| + |&n|) < oo. In fact, the 
general argument failed to capture Jacobi polynomials (whose clock es- 
timates were earlier obtained by Vertesi and Szabados [IHl 1121 UHl IM] ) 
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where separate arguments are needed. In |42ij, it was decided not to 
deal with asymptotically periodic OPRL since the derivative arguments 
looked to be tedious. 

Note. The papers [HI |12l Ej discussed asymptotics of OPUC when 
the Verblunsky coefficients decay exponentially. We have discovered a 
paper of Pan ^37] whose results overlap those in j5T1021E3- 

The key realization of this paper is that there is a more efficient 
way to eliminate pathologies like those in (|1.14j) . Namely, we will seek 
a priori lower bounds on eigenvalue spacings. If we find any 0{l/n) 
lower bound, that implies the rescaled p„ of (jl.lHj) has at most one 
zero near any zero of /oo- Such lower bounds are not new in suitable 
situations. Erdos-Turan already have such bounds if the measure is 
purely absolutely continuous in an interval with a.c. weights bounded 
away from zero and infinity. These ideas were developed by Nevai j^l] 
and Golinskii |T3]. Under suitable hypotheses on the transfer matrix, 
lower bounds are known in the Schrodinger operator community; see, 
for example, Jitomirskaya et al. [21] • {Note: They use "equal spacing" 
for 0(l/n) lower bounds and do not mean clock behavior by this term.) 

While we could have used these existing bounds in the proofs of The- 
orems 11.31 and 11.41 below, we have found a new approach which allows 
us to also prove Theorem 11.51 below, and this approach is discussed in 
Section 121 With these lower bound ideas, we can prove: 

Theorem 1.3 (= Theorem 13. 2j) . Suppose that 

oo 

^(|a„-l| + |6„|) <cx) (1.15) 

n=l 

Then there is uniform clock behavior on each interval [—2 + £,2 — e] 
for any e > 0. 

Theorem 1.4 (= Theorem 13. 4j) . Suppose that 

lim a„ = 1 lim &„ = (1-16) 

n— >oo n— >oo 

oo 

^(|a„+i - a„| + - 6„|) < OO (1.17) 

ra=l 

Then there is uniform clock behavior on each interval [—2 + e, 2 — e] 
for any e > 0. 

Remarks. 1. Theorem 11.41 implies Theorem 11.31 but we state them as 
separate theorems since the proof of the ffist is easier. 
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2. We will also prove results of this genre for perturbations of peri- 
odic recursion coefficients and for where n + 1 in (jl.lTj) is replaced by 
n + p for some p. 

3. We also obtain results (see Theorem 13 .Sp near ±2 if 

oo 

J2HK - l\ + \bn\) <oo (1.18) 

n=l 

Using our strong lower bound, we will also prove the following: 

Theorem 1.5. Let al^\b^^ be ergodic Jacobi parameters. Let Eq be 
such that 

(i) The Lyapunov exponent 7(-E'o) = 0. 

(ii) The symmetric derivative of v exists at Eq and is finite and 
nonzero. 

Then there exists C > so that with probability 1, 

hminf n[z«(Eo) - zI'\E,)] > C (1.19) 

n— >oo 

This result is especially interesting because it implies that the zeros 
cannot have Poisson behavior. It implies that Poisson behavior and 
ergodicity require 7(-Eo) > 0. We will say more about these issues in 
Sections HI and [TTJ 

While our initial interest in a priori lower bounds came from clock 
theorems and we could have finished the paper with Section |3] if our 
sole purpose was to prove Theorems ll.3Hl.5l it seemed natural to also 
consider upper bounds. Moreover, in looking over the upper/lower 
bound results in the OP literature, we realized one could get more 
from these methods, so we discuss that also. 

Broadly speaking, we have two sets of results and methods. The 
methods rely on either transfer matrices with hypotheses on recursion 
coefficients or on OP methods with hypotheses on the measure. We 
believe that the results are of interest to both the Schrodinger oper- 
ator and OP communities. Because the OP methods are unfamiliar 
to many Schrodinger operator theorists and are easy to prove (albeit 
very powerful), we have included an appendix with some major OP 
methods. 

The detailed plan of the paper is as follows: In Section |2l we prove 
our a priori lower bounds involving the transfer matrix (or, more pre- 
cisely, the growth of subordinate and nonsubordinate solutions). In 
Sectional we prove a variety of clock theorems, including Theorems 1 1.31 
and 11.41 In Section HI using ideas of Deift-Simon ISj , we prove Theo- 
rem 11.51 In Section 13 we obtain upper bounds on eigenvalue spacing 
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using the transfer matrix. Section IHl discusses using suitable Priifer 
angles to control spacing of zeros. 

In Section [71 we begin our discussion of OP methods with a technical 
result on bounds on 'w~^ and bounds of the Christoffel function. 
These bounds, which we will need for examples later are local versions 
of some bounds of Geronimus ^1] with a rather different method of 
proof. In Section we discuss upper bounds on eigenvalue spacing 
using OP methods and, in particular, find a remarkable lower bound 
on the density of states that is a kind of microlocal version of some 
bounds of Deift-Simon 0. In Sectional we discuss lower bounds on 
eigenvalue spacing. The methods in Sections |H1 and IHl are borrowed 
from Erdos-Turan [HIO], Nevai [Mj, and Golinskii ^1], but we show 
how to localize them and squeeze out stronger results. In Section 
we briefly discuss the analogs of our results for zeros of POPUC, and in 
Section [TT] discuss a number of examples, counterexamples, conjectures, 
and questions. 

It is a pleasure to thank L. Golinskii, S. Jitomirskaya, R. Killip, 
P. Nevai, and M. Stoiciu for useful discussions. This research was 
begun during B. S.'s stay as a Lady Davis Visiting Professor at The 
Hebrew University of Jerusalem. He would like to thank H. Parkas 
for the hospitality of the Einstein Institute of Mathematics at The 
Hebrew University, and E. B. Davies and A. N. Pressley for the hos- 
pitality of King's College, London where this was work was continued. 
Y. L. would like to thank G. A. Lorden and T. A. Tombrello for the 
hospitality of Caltech, where this work was completed. 

Percy Deift has long been a player in spectral theory and more re- 
cently, a champion for orthogonal polynomials. In particular, this pa- 
per exploits the work of Deift-Simon [HI. It is a pleasure to dedicate 
this paper to Percy. 



2. Variation of Parameters and Lower Bounds via 
Transfer Matrices 

Our goal here is to use variation of parameters to study eigenvalue 
spacing. Variation of parameters has an ancient history going back 
to Lagrange |15j and it was extensively used to study variation of so- 
lution with change in potential, for example, to study asymptotics in 
tunnelling problems fH]- The usefulness of the method as a tool in 
spectral theory goes back at least to the work of Gilbert-Pearson jl^ 
(also see f2ll2S|) with significant later contributions by Jitomirskaya- 
Last UHl dl 120] and Killip-Kiselev-Last It is essentially their 
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equation we will use although, interestingly enough, the earlier appli- 
cations are to E, E' with E eM. and E' = E + ie, while our application 
is to E, E' both in M. 

Given £^ G C, we consider solutions of 

dnUn+l + {K - E)Un + a„„iU„_i = (2.1) 

for n = 1,2,.... Here {a„,6„}^^ are the Jacobi parameters of the 
measure we are considering and 

ao = 1 (2.2) 
For 9 G [0, tt), we denote by Un{E, 6) the solution of (|2.1|) with 

Uq{E, 9) = sm{e) ui{E, 6) = cos(^) (2.3) 

In particular, 

Un{E,e = 0)=pn-iiE) n = l,2,... (2.4) 

The transfer matrix is defined by 

F^ _ f^n+i{E, 9 = 0) w„+i(E, 9 = 7r/2)\ 

^ " V «n(^, e = 0) u^{E, 9 = Ti/2) ) ^^-^^ 

so for any solution of ()2.H) . 

^"-')^T(n,E)M (2.6) 

Let K{n, m; E) be the kernel 

K{n, m- E) = m„,(E, 0)u„(E, 7r/2) - Un{E, ti/2)u.^{E, 0) (2.7) 
Define the operator Al{E) on C-^ = {f. = {f/tj^^i | tifc G C} by 

n 

(Ai(E)w)„ = K{n, m- E)vra (2.8) 

m=l 

(note that K{n,n) = 0, so the sum also goes to ri — 1). The following 
summarizes results from [THl UHl 1201 121| : 

Theorem 2.1. Let Wn{E),Wn{E') solve (|2.ip forE,E', and suppose 
WoiE) = woiE') wiiE) = w,{E') (2.9) 

Then 

n 

Wn{E') = Wn{E) + {E' -E)J2 K{n, m; E)w^{E') (2.10) 

m=l 

that is, 

w.{E') = w.{E) + {E' - E){Al{E)w{E')). (2.11) 
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Moreover, with || ■ ||hs the Hilhert- Schmidt norm on and defined 
by 

/ L X 1/2 

ll^-IU= El^^-I' (2.12) 

^ n=l ^ 

and ( ■ , ■)l the associated inner product, we have that 

PL(i?)||^s=lk(^,0)||i||«(E,7r/2/)||i-|KE,0),«(E,7r/2)),p 

(2.13) 

= max \\u{E,e)\\lmm\\u{E,e)\\l (2.14) 

u a 

In particular, 

L-l 

PL(i?)||<$^||T(j,E)f (2.15) 
i=o 

Remarks. 1. That (j2.10|) / ()2.11|) hold is either a direct calculation ver- 
ifying the formula or a calculation obtained by expanding {^^^^j^-^ ^) in 
terms of ("^^(ijf ) for ^ = 0, 7r/2. ^TH^ is a direct calculation from 
()2.7p and ()2.14|) is a clever observation 24 . 

2. Clearly, ||T(n,E)Qf = \u„.+im^ + > Km', 
so ||M(E,0)||i < RHS of ^A^. Similarly for ||m(E, 7r/2)||, so 
\\u{E,0)\\L\\u{E,7r/2)\\ < RHS of fT^ . so (jTTTH) implies fT^ . 

3. ||T(r?,) II measures the growth of the fastest growing solution, so the 
RHS of (ITTH|) in fact measures (max^ \\u{E,9)\\l) and thus, by TTPSi . 
one could place min^ ||m(-E', 6')1|/ maxe ||m(-E', 6')|| in front of the RHS of 

Here is the key lower bound: 
Theorem 2.2. LetE',E" be two distinct zeros ofpL{x). Then 

\E' - Eo\ + \E" - Eo\ > \\Al{Eo)\K (2-16) 

In particular, 

\zi'-'\Eo) ~ z[-'\Eo)\ > \\A,m)\\-' (2.17) 

Proof. By dTTTl) . for Ej = E', E", 

p.-iiE,) = p..im) + (E, - Eo)ALiEo)p.-iiE,) (2.18) 

To say Pl{E') = pl{E") = says E,E" are eigenvalues of Jl-f and 
p._i(-) are the eigenvectors. So, by orthogonality of eigenvectors, 

{p..^iE'),p._,iE"))L = (2.19) 
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By interchanging E' and if necessary, suppose 

\\p.-^{E')\\L>\\p.-r{E")\\L (2.20) 

Take ()2.18|) for E' and E" and take the inner product with p._i{E') 
and subtract to get 

\\p.-i{E')r < \E'-Eo\ \{p.,^{E'),AL{Eo)p.-^iE'))\ 

+ \E" - Eo\\{p.^^{E'),AL{Eo)p.-^{E"))\ 

<{\E'- Eo\ + \E" - Eo\)\\Al{Eo)\\ \\p.-i{E')r 

on account of ()2.2()|1 . ()2.1(')|) is immediate. 

f^n^ follows from TIU^ and z^'^) < < □ 

While our main applications are to clock theorems and Poisson sta- 
tistics, there is a universal tunnelling bound. 

Theorem 2.3. Let J^.p he a finite Jacobi matrix with a_ = inf 
a+ = maxa„, (3 = max6„ — min6„. Let 

J = aZ'liP + 2a+)2 + al + ll^/^ (2.21) 

Then any pair of eigenvalues, E, E' , of Jl-f obeys 

7^-1 

Remark. This bound is exponential, ~ 7~^^, for L large. 



E-E'\> (2.22) 



Proof. Adding a constant to 6„ does not change eigenvalue differences, 
so we can suppose that 

max6„ = — min6„ = ^ (2.23) 

Then any Eq in the convex hull of spec( J^^i?) obeys |-Eo| ^ f + 2«+, so 
|-E'o"~^ri| ^ /3 + 2a_|_. Thus, 7 is an upper bound on the Hilbert-Schmidt 
norm of 

1 (E - hn+i -an 
an+1 V 1 
and so on the norm. 

Pick E between E' and E". It follows that 

\\T{j,E)\\<j^ 

so ^ni^ follows from (jTTTI) and (IZTH|l . □ 

One can also use our proof to see that one cannot have too many 
zeros near E^. 
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Theorem 2.4. Define 5l,„(£'o) to be the smallest number so that 

i^{zf I \zf -E^\<5} <n-l (2.24) 

Then, for n > 2, 

h,n{Eo) >\\\Al{E^)\\^^^ (2.25) 

Remarks. 1. If one has strong clock behavior, „ ^ cnjL for n fixed 
and L large, so \Jn is worse than one expects in nice cases. 

2. Our proof shows that ()2.17|) can be "improved," if || ■ ||hs is used, 

to 

- z^,:'\e,)\ > \\A,iEo)\\^l (2.26) 

Proof. There are at least n zeros, zi, . . . , Zn, itl {z \ \z—Eq\ < 6L,n}- Or- 
der them so that if (pj^m = Pm{zj), then \\ipj\\L > ilv52|U > ■ > WfuWi- 
Let (fj = (fj/\\(fj\\L. Then the argument in the proof of Theorem 12.21 
says that, for j < k, 

l<[hA\An\ + \A,k\W 
<26lJ\A,,\' + \A,,\') (2.27) 

where 

Ajk = {Vj,AL{Eo)^k) (2.28) 
and ()2.27|) comes from (a + 6)^ < 2(a^ + 6^). Summing over all pairs 
and noting that l^np occurs n — 1 times, we find that 

n(n — 1) ^r. , X ,, , ,,9 

which is (jT^ . □ 

3. Clock Theorems for Bounded Variation Perturbations 
OF Free and Periodic OPRL 

The basic result from which we will derive all our clock theorems pre- 
supposes the existence of a complex solution to the difference equation 
1)2.11) for which we have precise information on the phase. The model 
is the Jost solution which is asymptotic to e*^*^-^)", where E = 2cos(^), 
< ^ < vr. 

Theorem 3.1. Let {an,bn}'^^i be a set of Jacobi parameters and A a 
closed interval in R. Suppose there exists a solution Un{E) of ()2.H) for 
e A which obeys 

(i) 

uo{E) > lmui{E) > (3.1) 
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(ii) 

Un{E) = \un{E)\eM^[na{E) + (3.2) 



where a is with 



(icy 

_^>0 allEeA (3.3) 
dE 



each (3n continuous on A, and 



lim 

n— »oo 



sup |/3„(E')-/3n(^)| 
\E-E'\<l/n 



(3.4) 



(iii) The transfer matrix T„{E) of (j2.5|) obeys 



T= sup ||T„(E)|| < oo (3.5) 

n,E£A 



Then the density of states exists on A, 

du(E) = --^dE (3.6) 

TT dE 

and there is clock behavior uniformly in A. 

Remarks. 1. a is, of course, a rotation number and ()3.6|) an expression 
of the connection between the density of states and rotation numbers; 
see Johnson-Moser ^22j. 

2. ()3.4|1 imphes f3n is irrelevant for eigenvalue spacing comparable 
to 1/n. To control possible spacings with AE small compared to 1/n, 
one needs some Lipschitz control of that is, 

.up lfl.(^')-ftf)Uo 

Ej^E' n\E - E'\ 

which is where differentiability of /?„ and so moment conditions on 
{an,bn} came into [41J. We avoid this by using ()3.5p to get a priori 
bounds. 

3. (j3.4j) implies the same if 1/n is replaced hj A/n for any fixed A. 
Define UA) by 

UA)= sup \(3n{E') - f3niE)\ (3.7) 

E,E'gA 
\E-E'\<A/n 

4. ()3.4p is implied by an equicontinuity assumption, for example, 
uniform convergence of to a continuous limit. 

Proof By fj3.1|) . m„ and independent solutions of ()2.1|) and so 

cannot vanish at any points. Moreover, 



Pn-i{E) = A{E)un{E) + A{E) Un{E) (3.8) 
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where 



A{E) 



(3.9) 



W{u,u) 

Here, W is the Wronskian. Given sequences, fn,gn, 

Wnif, g) = anifn+l9n " fn9n+l) (3.10) 

is constant (call it W{f,g)) if f,g both solve ()2.1|1 . Since p_i = 0, 
Po = 1, gq = 1, we have W{p.^i,u) = Uo{E), and clearly, W{u,u) = 
22 Im('Ui'Uo) = 2iuQlmui, so 



A{E) 



1 



(3.11) 



2 IniMi 

is pure imaginary. Thus, ()3.8p says Pn-i{E) vanishes if and only if u 
is real, that is, by ()3.2p . 

p„_i(E) = ^ na{E) + l3n{E) = kn 

Let 



{k e z) 



11 = mm 



da 
"dE 



(3.12) 
(3.13) 



Pick A'o so that n > Nq implies 

1] 



>4Cn 



(3.14) 



This can be done since Cn{A) as n ^ cxo by hypothesis. Since 
Cn{A.) is increasing in A and Cn{A + B) < Cn(^) + Cn(-B), we have 
Cn(a;/3) < U[xm + Ux(3-[x]f3) < {[x] + l)U(3) < 2xU(3), iix>l, 
(jmil implies 



q> m> 2Cn(g) 



(3.15) 



This in turn implies 
1 



E-E'>^^ [n[a{E') - a{E)] + [/3„(i?') - f3^{E)]] > - i]\E' - E\ 



llT^ 



By fl2.17|l and ()2.15|1 . any two successive zeros obey 

1 



\E' -E\ > 



llT^ 



(3.16) 



(3.17) 



Thus, fl3.17|) implies that for n > Nq, any two solutions of ()3.12|) have 
distinct values of k. We also see from ()3.16|) and continuity that if E 
is a solution of ()3.12|) . there is another solution in {E — E) and it 
has the next larger value of k (i.e., k + 1). 
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Subtracting ()H.12j) for two successive values of ()H.12j) and using 

Cn(2vr/?7) — > 0, we see that 

sup \n{a{E') ~ a{E)) - n\ ^ (3.18) 

E' <E successive 
eigenvalues in A 

Given the uniformity of convergence of the difference quotient to the 
derivative, ()3.18p imphes that 



sup 

E'<E successive 
eigenvalues in A 



(3.19) 



This imphes the density of states exists and is given by ()3.6p and 
that one has uniform clock behavior. □ 

Theorem 3.2 (= Theorem II. 3p . Let {a„,&n}5JLo ■set of Jacobi 
parameters obeying 

oo 

Y,\bn\ + \an-l\<oo (3.20) 

n=l 

SO esssupp((i/i) = [—2,2]. For any e > 0, we have uniform clock 
behavior on [—2 + e,2 — e]. 

Remarks. 1. This includes Jacobi polynomials (rescaled to [—2,2]) for 
which 

\bn\ + \an - 1| = 0(n"^) 
2. Of course, the density of states is the free one. 

Proof. It is well-known (see, e.g., j2Slini) that when (|3.2(jp holds, there 
exists, for all x G (—2,2), a solution Un{x) so that if 

z + z-^ = x (3.21) 

(i.e., z = e*^ with x = 2cos0 and 6 G (0,7r)), then 

z-'^Unix) 1 (3.22) 

uniformly on compact subsets of (—2, 2). Moreover, Un{x) is continuous 
on (—2, 2) for each fixed n. By evaluating the Wronskian near n = oo, 
we see 

WiUn,iLn) = Z- (3.23) 

Thus, if 

Un{E) = ^^ (3.24) 
uo{E) 

then dUl) and ^ hold. If 

uo{E) = |uo(^)|e'^°°(^) (3.25) 
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then we have that 

Pn{E) ^ P^{E) (3.26) 

uniformly on [—2 + e,2 — e]. By continuity, ()3.4|1 holds. 
a is given by 

a(E) = arccosf ^ ) (3.27) 



^ 2 
so 

1 da 1 1 



(3.28) 



IT dE TT - 

and dsn holds. 

Finally, standard variation of parameters about z^, z^" shows that 

holds for each A=[2 + e,2-e]. 
Thus, Theorem 13.11 applies, and we have clock behavior. □ 

In the above, we used the fact that /3„ — > P<x uniformly to obtain 
fl3.4|) . In the bounded variation case, we will instead use: 

Lemma 3.3. If (3n = ^n'^ + where (3^^^ is and 

"^^^ (3.29) 



1 

— sup 



dE 

('2') (2,') JT^TTT 

and l3n — > (3oo uniformly, then ()3.4j) holds. 
Proof. Immediate since E,E' G A implies 



dx 



\Pr.{E) - Pn{E')\ < \E~E'\ sup 

EGA 

Theorem 3.4 (= Theorem II. 4|) . Let {an,&n}^o ^ J(^cobi 
parameters with 

a„ ^ 1 &n (3.30) 

and 

oo 

^(|a„+i - a„| + - 6„|) < oo (3.31) 

n=l 

then for any e > we have uniform clock behavior in [—2 + £, 2 + e] . 

Remark. Again, the density of states is the free one by ()3.30p . 

In order to prove this theorem, we need the following result: 

Theorem 3.5. Let Bn{9) depend continuously on 6 E I, a compact 
subinterval of (0, tt), and suppose 

oo 

sup Y \\Bn+i{e) - Bn{e)\\ < OO (3.32) 
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and 

Bn{e) ^ B^{e) (3.33) 

uniformly where -Boo(^) has eigenvalue e^*^. Explicitly for 1^00(6*) con- 
tinuous and invertible: 



Bue) = v^{e)[^^ ^:..jKo(^)-^ (3.34) 

Suppose each Bn{9) has eigenvalues e"^^"^"^^^ with ^Pnid) G (0,7r). Let 

fnie) = Bn{e)...B,{e) (3.35) 

Then 
(i) 

sup \\fn{e)\\ < 00 (3.36) 
eei,n 

(ii) There exists Soo{d) depending continuously on 9 so that 

\\fn{e) - v^{e)Dn{e)s^{e)\\ ^ o (3.37) 

where 

Dn{0) = i^' Q ^_ '„(,)J (3.38) 

with 

n 

r]4e) = Y,M0) (3.39) 

v=i 

Remark. See the notes to Section 2.1 of [IHl for a history of resuhs on 
bounded variation. 

Proof. This is a strong version of Kooman's theorem |29j. In Sec- 
tion 12.1 of 40j, (j3.36p is proven, and in the notes to that section, it 
is noted that ()3.37|1 holds. In those notes, there is no Voo because the 
analog of D„ is not diagonal but can be diagonalized in a basis where 
Boo is diagonal. □ 

Proof of Theorem Y^. Define 

^ f (a„a„+i)^^/^(2cos(e) - 6„+i) -(a„/a„+i)^/^A ^^^^ 

SO the transfer matrix at = 2 cos 6 is 

Ue) = (^) "\^{e) l"-^) B,{6) 

= a-li'Ue) (3.41) 
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Since det(i?„) = 1, Bn{6) has eigenvalues e^*'''", ipn G (0,vr), if and 
only if 

(a„a„+i)-^/2(2cos0 - = 2cos(v^„) G (-2,2) (3.42) 

and this holds uniformly for 9 E (5, tt — £:) and n > Nq for some fixed 
Nq. Thus, for n > Nq, we define 

n 

vr.{e) = (3-43) 

j=No 



and define 



and 



Dnie) = ( - Q ) (3.44) 



'S'oo(^) = Sooid; Bn - ■ ■ BN^^)B]y^_i . . .Bi (3.45) 

where Sqo{6] Bn ■ ■ ■ Bnq) is the Soo for the sequence B^^, Bmo+i, ■ ■ ■ ■ 
Thus, fP?37|l and (HOTIl show that 

- a-l{'V^{e)Dnie)SUe)\\ ^ O (3.46) 

We now proceed to construct a solution u obeying the hypothesis 
of Theorem O Pick a two-vector x{e) by x{e) = e*^"(^)^oo(^)"^ (J) 
where ipo is chosen below and m„ by 

T„(^^)x(^^) = f'^-^;!^^! (3.47) 



then (j3.46|) says that 

un{e) = \un{e)\exp{i[7]n{e) + pl^\e)]) (3.48) 

where 

P^:\0)^Pg\9) (3.49) 

uniformly. 

(2) 

Here /?^^(^) is v5o(^) plus the phase of the 21 element of ^0(6*), 
and so i{r]n{9) + (3oo {d)) is the phase of the lower component of 
Voo{d)Dn{0)Soo{())x{9). Since is not real, mo(6') 7^ 0, and so 
V9o(6') can be chosen so that ^0(6*) > 0. 

Since ^jiO) — > 6, rjn+i — rjn ^ 0, and thus the imaginary part of the 
Wronskian of u and u is positive, so u obeys (j3.1|) . 

By 

dipn sin(g) ^1/2 c-n^ 

ai? sm.[ipn) <jE 
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SO 

dE dE ^ Vsin(<^n) / 



(3.51) 



j=No 

Thus, if 

de 



/3« = r7„-n^ (3.52) 



then 



1 / N,l^ [ sin(g) /2 \\de 

- = ^ ~ — V (a„a„+i) / - 1 ^ (3-53) 

n oE \ n n \m\\}Pn) ) ) dE 

converges uniformly to zero, since 

1 (a„a„+i)-^/' 1 (3.54) 



Sin(^) , T „ ^-l/2 



sin(v9„) 
uniformly in 6. 

Lemma 13.31 applies, so condition (ii) of Theorem 13.11 holds with 
a{E) = 9{E) = arccos(-|). We thus have clock behavior with den- 
sity of states the free one, that is, given by ()3.28|) . □ 

Now consider the periodic case, that is, 

(^n+p — (^n bn+p = bp (3.55) 

The spectrum now has bands (see, e.g., |20|). For any E in the interior 
of the bands, there is a Floquet solution with Un+p = e^'^^^^w^ with 
7(E) e (0,7r) and || < 0. Thus 

ULp+r = |w.|e^(^->(^)+^'°°') (3.56) 

where (3r°°^ is the phase of Ur- Theorem 13.11 applies with a = 'y{E)/p 
and 

f3Lp+r{E)=Pi^\E)-^-^{E) (3.57) 

There are only r such functions so ()3.4j) holds, and we recover the zero 
spacing part of Theorem 2.6 of ^1^. 

If an \ bn^ are periodic and a„ = + dan, bn = ^ + 5bn and 

oo 

^(|56J + |5a„|) < oo (3.58) 

n=l 

then one can construct Jost solutions on the interiors of the bands. All 
that changes is that ()3.57|1 is replaced by 

lim l3Lp+r{E) = f3i°^\E) - - -f{E) (3.59) 

L^oo p 
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SO Theorem Kill still applies. Similarly applying the ideas in the proof 
of Theorem 13. 41 we obtain a bounded variation result. Since it includes 
the ()3.58|) result, we summarize in a single theorem: 

Theorem 3.6. Let an\hn^ obey 

'^n+p — '^i ^ ^n+p — ^ (3.60) 

for some p. Let an, bn obey 

lim |a„-4°)| + |6„-6(,°)| = (3.61) 



n^oo 
oo 



^ (kn+p - a-nl + \bn+p - &n|) < OO (3.62) 
n=l 

Then, for any closed interval which is a subset of the interior of the 
bands {see Remark), we have uniform clock behavior with density of 
states —- 

p dE 

Remark. There are p closed bands, . . . , Bp, generically disjoint but 
perhaps touching in a single point (closed gap). By the interior of the 
bands, we mean Vf-^^BY^ which will be smaller than (U^^^Sj)™* if some 
gap is closed, that is, we must remove all of the gaps, including those 
that degenerate to single points. 

p need not be the minimal period, so we have that 

Corollary 3.7. Suppose 

a„ ^ 1 ^ 

and, for some p, ()3.62j) holds. Then, in any closed interval in {E = 
2 cos{6) \ p6 ^ 0, ^, • • • , ^~p^^ }, we have uniform clock behavior. 

As a final topic, we want to discuss zeros very near E = 2 when 

oo 

5^n(|a„- 1| + |6„|) < oo (3.63) 

n=l 

It should be possible to extend this argument to get uniform clock 
behavior in [—2,2], with a suitable modification to take into account 
the behavior exactly at ±2. When (j3.63j) holds, the Jost function, u, 
can de defined on [—2,2]; see, for example, the appendix to 0. If 
u{2) = 0, we say there is a resonance at 2, and if u{2) ^ 0, we say that 
2 is nonresonant. 

Theorem 3.8. Let {a„, be a set of Jacobi parameters obeying 

(jSSni). Define < ^J"^ < ^^"^ < ■ ■ ■ so E^^ = 2cos(^f ^) are the zeros 
of Pn{x) nearest to x = 2 and below. Then 
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(a) If 2 is a resonance, then 

nef - (j - i)7r (3.64) 

(b) If 2 is nonresonant, then 

nOf^ ^ JTT (3.65) 

Remarks. 1. The two simplest examples are the nomesonant a„ = 1, 
6„ = where 

p„(2cos^) = c„ (3.66) 
sm 6 

and the resonant a„ = 1 (^ > 2), Oi = -\/2, and 6n = 0, where 

p„(2cos6') = dnCos{n9) (3.67) 

2. Notice that (for simplicity, consider the nonresonant case) 

(„) („) (2j + l)7r^ 
-E,^i 

the eigenvalue space is 0{n~^) and not equal in E but has clock spacing 
in e. 

Remark. The key fact that at a zero energy resonance, the scattering 
phase is 7r/2 (mod vr) and otherwise it is (mod tt) is well-known in 
the continuum case, for which there is extensive physics literature; see, 
for example, Newton pUj . 

Proof. By the theorems found in the appendix to [HI (which codifies 
well-known results), when fl3.63|) holds, one has (0 eqn. (A. 27)]) 

b„(e*')| <C(n + l) (3.68) 

and the existence of a solution Un{e^^) with 

e-*"^u„(e^^) 1 (3.69) 

uniformly on dVi. Uo is called the Jost function and 

Wip._,,u)=Uo (3.70) 

We want to use ()2.14|) where there is a collision of notation, so we 
let v{9,Lp) be the solution at E = 2cos^ and boundary condition ip. 
Then fl2.14|) becomes 

Pl(2cos^)|| <max \\v{e,<^)\\Lmm \\v{e,<p)\\L (3.71) 

ip LP 

If uq{6 = 0) 7^ 0, we get one solution for 6 small, v{6,ipo{6)) which 
is uniformly bounded in 6 and n, and another solution v{6,Lpi = 0) 
(= p.-i) bounded by Cn. It follows that 

min \\v{e,y^)\\L < CL^/^ max ||i;(^,<^)|| < CL^'^ (3.72) 
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If uo{0 = 0) = 0, we start at n = 1 since mi(6' = 0) 7^ and construct 
the bounded and linearly growing solution that way (in essence, the 
two solutions in this case are and q, where q is the second kind 
polynomial), so ()3.72p still holds. 

We conclude, using Theorem 12.21 that when ()3.63p holds, then 

sup \E -E'\>A^ (3.73) 



E,E' successive zeros of p„ 
E,E'e\2-e,2] 



71^ 



Define (p{0) by 

Mo(e'^) = |Mo(e'')|e'^(') (3.74) 

Lf can be defined by ^ 7^ since uq{9) is then nonzero. We can pick 9 
continuous on (0,e). We claim 

fO mod vr if 2 is nonresonant 
<^(0) = lim <^(^^) = <^ ^ (3.75) 
ej-O I ^ mod vr 11 2 is a resonance 

Postponing the proof of this for now, let us complete the proof of the 
theorem. 

By fj3.70|) and W{u, u) = z — (since u ~ 2;"'), we see that 



Pn{2cose) = (3.76) 

-l"°^^")"""^^")isin(n^ + /?„(^)) (3.77) 



2sin^ 
where 

(3^M ^ -ifie) (3.78) 

as n goes to infinity uniformly in a neighborhood of 6 = 0. Thus, zeros 
of Pn{2 cos9) are given as solutions of 

nO + (3n{e) = JIT (3.79) 

In the resonant case, since 



sup 

I I — n 



^ (3.80) 

mod TT 



there is at least one solution asymptotic with 

nef ^ (j - Dtt (3.81) 

If there were multiple solutions for some j, we would have two zeros 
with 0<9' <9, 

n{d-e')^o e <- 

n 
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for n large (vr can be any number strictly larger than 7r/2), 
|2cose'-2cos^| < |^'-e|sin(^) 



n J \n 

violating ()3.73|1 . Thus, there are unique solutions and ()3.64|1 holds. 
In the nonresonant case, ()3.79|1 holds, but instead 

sup \(3n{0)Uod.^0 (3.82) 

which proves existence of solutions with 

nef ^ jTT (3.83) 

for j = 1,2, ... . We must prove uniqueness for j > 1 and nonexistence 
for j = 0. 

The uniqueness argument for j > 1 is the same as in the resonant 
case. To show no solution with j = 0, we suppose that J has m 
eigenvalues above E = 2 (by Bargmann's bound 17^, the number is 
finite). Let J(A) be the Jacobi matrix with 



a„(A) 

bn + \ n < m + 1 
bn n > m + 1 



bnW 



It is easy to see that as A — oo, J(A) has at least m + 1 eigenvalues. 
So pick Ao = inf{A | J(A) has m + 1 eigenvalues in (2, oo)}. Then 
Ao > and J(Ao) has a resonance at 2. By the analysis of the resonant 
case, pii'°\x) has m zeros in (2, oo) and its (m + l)st zero asymptotic 
to 2 - (i If, which means Pn{x), whose zeros are less than those of 

pi^"'', cannot have a zero asymptotic to 6'^^'^ 0. 

That proves the result subject to ()3.75j) . In the nonresonant case, 
Mo(e*^) is continuous and nonvanishing at 6^ = 0, and mo(1) is real, 
(/}(0) = mod vr, so continuity proves the top half of p.75|) . 

In the resonant case, we note that a Wronskian calculation (see [HI 
eqn. (A. 49)]) shows that 

Im(Mi(e^^)Mo(e*^)) = sin0 (3.84) 

Since mo(1) = 0, mi(1) 7^ 0, and mi(1) is real, so 

liinlmf^^^^ ^0 (3.85) 



eio V sin 9 
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On the other hand, Mo(e~*^) = 'Uo(e*^), so Lp{0) + '^{—0) = mod vr, 
which means that any hmit point of (p{0) is a multiple of tt/2. This is 
only consistent with ()3.85j) if the limit is congruent to 7r/2 mod tt. □ 

4. Lower Bounds in the Ergodic Case 
Our main goal in this section is to prove the following: 

Theorem 4.1 (= Theorem II .Sj) . Letal^\b^^ be ergodic Jacobi param- 
eters. Let Eq be such that 

(i) The Lyapunov exponent '^{Eq) = 0. 

(ii) The symmetric derivative of v exists at Eq and is finite and 
nonzero. 

Then there exists C > so that with probability 1, 

hminf n[z«(Eo) - ^^^H^o)] > C (4.1) 

n—foo 

This is particularly interesting because of the connection to Poisson 
behavior, where: 

Definition. We say a probabilistic family of Jacobi matrices has Pois- 
son behavior at Eq if and only if for some A (normally A = density of 
zeros) we have that for any ai < Pi < a2 < /32 < ■ ■ ■ < Pi and any 
ki,k2,...,h e {0,1,2,...}, 



Prob(#<'z(^)(Eo)G 



n n 



km for m 



kr, 



n 

m=l 

(4.2) 



e 



Remark. Poisson behavior was proven in certain random Schrodinger 
operators by Molchanov and for random Jacobi matrices by Mi- 
nami [32]. See Stoiciu 071 EHl and Davies-Simon V for related work 
on OPUC. 

Corollary 4.2. Let a^n \ bl^^ be ergodic Jacobi parameters and Eq 

so that the symmetric derivative ofv exists at Eq and is finite. Suppose 

there is Poisson behavior at Eq. Then 7(-E'o) > 0. 

Remarks. 1. Basically, ()4.1|) is a rigid level repulsion inconsistent with 
Poisson behavior. 

2. Ergodicity is critical here. Killip-Stoiciu [27] have examples which 
are not ergodic for which there is Poisson behavior with 7 = 0. 
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Proof of Corollary Suppose first 7(-Eo) = so Theorem 14.11 ap- 
plies. For each ra, let fn{^) be the characteristic function of {uj \ 



^'\Eo) - zi'''\Eo)\ < iCn-i}. By dH)), /„(cu) ^ for a.e. u, 

fniuj)duj^O (4.3) 



n 



SO 



as 72 — oo. 

Clearly, if there is one z^^^ in [Eq — jCn~^, Eq] and one in [Eq, Eq + 
jCn~^], then /„(cij) = 1. Thus, by the assumption of Poisson behavior, 

hm j f.{u:)du> (^e-^^/^y 

which contradicts (j4.3j) . □ 

Our proof of Theorem 14 . II will use the complex solutions constructed 
by Deift-Simon [S] and the estimate Theorem 12.21 It is thus impor- 
tant to be able to estimate A in terms of any pair of solutions with 
Wronskian 1. 

Lemma 4.3. Letu^^\u^'^^ he any pair of solutions of (j2.H) with Wron- 
skian 1. Then the kernel K of (j2.7j) has the form 

K{n,m,E)=u^^\t^-u^^\^^^ (4.4) 

In particular, 

Proof. Noting that Un^Um — UmuiP = 0, we see that K is invariant 
under linear changes of the m's of determinant 1. This proves ()4.4|) . 
()4.5p follows as in Theorem 12.11 □ 

We need the following result of Deift-Simon jH]: 

Theorem 4.4. Let a^\b^^ be ergodic Jacobi parameters and let (i)- 
(ii) of Theorem \4 .l\ hold for Eq. Then for a.e. uj, there exists a complex- 
valued solution u{- ,uj) of 1)2.11) so that 

(i) The Wronskian of u and u is —2i. 

(ii) We have that 

n— 1 ^ 

limsup \n\~^ |M(j, w)\^ < 27r - — (4-6) 



Remarks. 1. In fact, \u{j,w)\ = \u{0,T^w)\ and E(|m(0,w)P) < oo, so 
by the Birkhoff ergodic theorem, the lim sup on the left of ()4.6|) can be 
replaced by a limit which is a.e. constant with a constant bounded by 
the right side of ()4.6|) 
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2. jH] states results for a.e. E with ■j{E) = 0, but the proof shows 
that what is needed is (i)-(ii). 

Proof of Theorem \4.1\ By ()4.5|) . taking u^^'> = u and u^"^^ = {2i)^^u, we 
see 

hmsupn-lA„(E)|| <7r^ (4.7) 



by ()4.(i|l . Thus, by Theorem 12. 2^ 



hminf n[z«(Eo)-4"'n^o)] > - 

n— >oo TT 



_dE_ 

which is dHH). □ 

5. Upper Bounds via Transfer Matrices 

Our goal in this section is to prove: 

Theorem 5.1. Let {an, &n}5JLi be a set of Jacobi parameters. For any 
bounded interval / C M, we have 

,^ [8esup^,,||T„(E)||](n;=,a,)^/" , ^, 
sup \E - E \ < (5.1) 

E,E' successive ^ 
zeros of pn in I 

Corollary 5.2. Let {an, bn}'^=i be a set of Jacobi parameters and I = 
[a,j3] a closed interval. Let 

y4 = supa„<oo T= sup ||T„(i?)|| < cxd (5.2) 

n n\E(iI 

Let c = 8eAT. Then for any 6 > 0, there is Nq so that if n > Nq and 
-Eg [a + 6, f3 — S] is a zero ofpn, then there are at least two additional 
zeros in \E — -, E + -], one above E and one below. 

Proof of Corollary 5.2. It is known that any point in spec(J) is a limit 
point of zeros, so for large enough Nq, there are zeros in [a, a + 6) and 
(/3 — S, P]. Thus, ()5.H) implies the result. □ 

Example 5.3. Let a„ = 1, 6„ = so 

sm((n + 1)^) 
p„ 2cosg = ^\ , ^ ^ 5.3 
smt^ 

and for n odd, Pn{0) = 0. The next nearest zero is at 6* = | — so 
at E ~ ^ (^(2cos^)|^^^ = -2). In this case, p„ at E = (0 = f ) 

is (1, 0, -1, 0, 1, 0, . . . ) for n = 0, 1, 2, . . . and g„ = (0, -1, 0, 1, . . . ), so 
T„odd = {^i V) so ll^n(0)l| = 1 and ||T„(E)|| ~ 1 for E near 0. Thus, 
the correct answer for the spacing is 27r ~ 6.3 and our upper bound is 
8e ~ 21.7, a factor of about 3.5 too large. □ 
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To get Theorem 15. H we will use 

Theorem 5.4. Let Q be a polynomial with all its zeros real. Let 
Q{Eq) = 0, Q'{Ei) = 0, with Eq < Ei andQ nonvanishing on {Eq, Ei). 
Then 

Proof. Since ^ log(x — Xq) = — |^_^p|2 < 0, we see that g{E) = 
\og{\Q{E)\/\E - Eo\) is concave on [Eq, Ei]. Note that g{Eo) = 
log\Q'{Eo)\ and that the tangent to g at Ei is 



\Q{E: 



Thus, 

\Q'{Eo)\< 

which is dniH). □ 
The Q we will take to get (j5.1j) is not P, but the discriminant 

A„(E) = Tr(r„(E)) (5.6) 
associated to the periodic set of Jacobi parameters 

^mn+q ~ ^mn+q ~ ^q q = 1, . . . , Tl] m ^ (5.7) 

We have: 

Lemma 5.5. The zeros ofpn-i and An interlace. Thus, if Ei < E2 < 
E^ are three successive zeros of An, thenpn^i has two zeros, E and E' , 
with 

\E - E'\<\E^~ Ei\ (5.8) 



Proof. We need the analysis of A„ as a periodic discriminant 
A„ has n bands given by ai < /3i < 0^2 < < ■ ■ ■ < a„ < /3„ 
and bands [aj,[3j\. A^^{—2,2) = U^^^{aj, Pj) and, in particular, 
the zeros of A lie one per band. Pn-i has one zero in each gap 
[Pi, ^2], [/52, ^3], • • • , [Pn-1, «n]- That gives us the interlacing. ()5.8|) is 
an immediate consequence of this interlacing. □ 

To get a bound on A^ at its zeros, we need a bound on the rotation 
number for ergodic Schrodinger operators found by Deift-Simon 8J. 
This rotation number is 

a{E) = n{l-i^{-oo,E)) (5.9) 

where u is the density of states. Thus, a runs from vr to as ii^ runs 
from minspec(J) to maxspec(J). cos(q;) runs from —1 to 1. 
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Proposition 5.6. In the periodic case, on spec(J), 



dcos{a{E)) ^ 1 
dE - 2 



n 



-l/n 



(5.10) 



Remarks. 1. (1.2) of jH] is an integrated form of (j5.1(Jj) . We can take 
derivatives since spec (J) is a union of intervals. 

2. Deift-Simon assume aj = 1. By using the modification of the 
Thouless formula for general aj, it is easy to see their proof yields 

(inni). 

3. In the free case (a„ constant, 6„ = 0), one has equality in (jS.lUp . 
Proof of Theorem \5.1\ An{E) is related to a{E) by 

A„(E) = 2cos(na(E)) (5.11) 

Thus, at zeros of A,„ where cos{na{E)) = 0, we have \sm{na{E)) \ = 1. 
So at such points. 



dE 



2n^>2nsm{a{E))^ 



In 



> n 



dE 
dcos{a{E)) 
dE 

n -l/n 



by (ICTH) . 

Since A„ = Tr(T„), |A,„| < 2||T„||, so ()5.4j) becomes 



\Ei - Eo\ < 2e 



\Tn\m 



^l/n 



n 



(5.12) 



Between the zeros of A are two {Eq, Ei)-tjpe intervals and so, be- 
tween the first and third of three zeros are four such intervals. fl5.8|) 
and (jKT^ imply (jSH)). □ 

6. Prufer Angles and Bounds on Zero Spacing 

There are various possible Priifer angles. We will exploit one that is 
ideal for studying the energy dependence of zeros of 

Proposition 6.1. Fix Jacobi parameters {an,bn}'^=i- For each n = 
1,2,..., there is a unique continuous function 6n{E) determined by 

PniE) 



tan(^„(E)) 
;im en{E) 



Pn-l{E) 

~2 



(6.1) 
(6.2) 
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Moreover, 

dE a„(p„_i(E)2+p„(E)2) ^"-^^ 

Proof. pn{E) / pn-i{E) — oo as E ^ — oo and the ratio is continuous 
on M U {00} , so existence and uniqueness are immediate, as is ()6.2|) 
since tan(— 7r/2) = —00. Note next that since ^ arctan(?/) = jq^, we 
have 

dOn ^ Pn-lP'n " PnP'n-l 

dE pU{1+pI/pL,) 
SO that follows from the CD formula ()A.3j) . □ 

Remark. ()6.3p is closely related to well-known formulae relating the 
derivative of projective angles in SI[J(1, 1) or §L(2,M); see Theo- 
rem 10.4.11 of jini- These formulae have been used by Furstenberg, 
Carmona, and others; see, for example, Carmona [5 or Simon 



The transfer matrix, Tn{E), is a 2 x 2 matrix of determinant l/a„+i, 

so 

||T„(E)-i =a„+i||T„(E)|| (6.5) 
Thus, since {pj+ijPjY = Tj(_E')(l, 0)*, we have that 

ajUTAE)\r < P,+i{Ef+pAEf < WUEW (6.6) 

Since 

n—2 n—1 n—2 

iE(^l+^?.i)<E^?<E(^?+^?.i) (6-7) 

i=o j=o i=o 

(16. 3j) immediately implies 

Theorem 6.2. Let 

tn{E)= sup {I + a%,)\\T,{E)\\^ (6.8) 

0<i<n-l 

Then 

2«„(E)L.(g) -''"^-'""'^"'-'<^' (6.9) 

This, in turn, implies 1/n upper and lower bounds on zero spacings 
sufficient for what we needed in Section IHl 

Theorem 6.3. //A is an interval in M on which = sup^g^^ (1 + 
al^^)\\Tn{E)f then 

inf \E-E'\>^ (6.10) 

E,E' successive zeros of p„{E) 
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and 

sup |E-E'|<^^^ (6.11) 

E,E'€A n 
E,E' successive zeros of pn{E) 

Moreover, if A = [a, j3], pn has zeros in ^ /5] and [a, a+ ^""^'^^ ] 

once |p — q;| > — 

Proof. Since 6'„ is monotone in and Pn{E) = if and only if 6n{E) = 
in for some £ G Z, we have at successive zeros, E < E', that 

fE' in 

'n = en{E')-er,{E) = -^dE (6.12) 

()6.1()|1 / ()6.1H1 are then immediate from ()6.9p . The final assertion comes 
from the fact that 6{Ei) — 6{Eq) > it implies that tan(6'(i?)) has a zero 
m[Eo,E,]. □ 



7. Relations of the Weight to the Christoffel Function 

The previous sections were dominated by the transfer matrix. In 
this section, we shift to the weight where the CD kernel (see ()A.1|) ) 
will play a major role. This section is a technical interlude: a detailed 
result that will be useful in the analysis of examples in later sections. 
Our main result in this section is 

Theorem 7.1. Suppose that 

djji = w{x) dx + djis (7.1) 

with dfj,s singular, and that for some xq, a > 0, and some r > 0, 

rXQ+a 

/ wix^dxKoo (7.2) 

J XQ—a 

Then 

K4xo,xo)<Crn'+'-~" (7.3) 
where Cr only depends on r, a and the integral in ()7.2j) . 

This result generalizes one of Geronimus (see Remark 3.3 and 
Table II]) in two ways. His estimate is on |</?nP not Kn and, more 
importantly, his estimates require global estimates on w in the context 
of OPUC rather than just our local estimate. One reason we can go 
beyond Geronimus is that he uses the Szego function and we just use 
the Christoffel variational principle. Another reason is that we have a 
powerful result of Nevai j^l] : 
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Proposition 7.2 For any p in (0, oo), there is a constant Dp 

so 

J \7Tnix)\Pdx>Dpn-'\nn{0W (7.4) 
for any polynomial Tin of degree n. 

Proof. Since this is a special case of Nevai's result which depends on 
several arguments, for the reader's convenience, we extract exactly 
what is needed for ()7.4|) . 
Let 

dfioix) = - {1 - x^)-^/^ dx (7.5) 

TT 

on [—1, 1], so if X = cos(6'), then 

d9 

dfiQ = — 

TT 

on [0, vr] which implies, as is well known, that the OPs for ()7.5|1 are 
given by 

Pn{cos{e)) = r ^ = ° (7.6) 

I yz cos(n6') n > I 

the Chebyshev polynomial (of the first kind). 
It follows that for x = cos{6) G [—1, 1], 

n 

Kn{x, x; dfio) = 1 + ^2 cos^(jn) < 2n + 1 (7.7) 
i=i 

Thus, by (pTTj) . 

sup Kix)\' < {2n + 1) / Kix)\'il - x2)-V2 _ (7.8) 



hr <1 



TT 



If m is an integer, (vr^)™ is a polynomial of degree at most mn, so 
()7.8p implies 

/I r/r 
\7ln{x)\^"'{l-xY'^^— (7.9) 
-1 ^ 

If 2m — 2 < p < 2m, we write 

. 2m— p 



\x\<l 

to deduce 

sup |7r„(x)|^ < {2n 



\x\<l 



TTnl^"" < |7r„,|*'( sup |7r„(x)|) 

-1/2 



7r„(x)r(l-x2)-V2_ (7.10) 

TT 
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Given p, pick i from 1,2, .. . so ip > | and apply (j7.1U|) to the 
polynomial 7r„(x)(l — x'^Y which has degree n + 2i and get 



k„(0)r < sup \{l-xyTin{x)\^ < { 2{n+2(i) 

\x\<l 



+ 1 



TT 



since (1 — x'^) 



(7.11) 



-5 < 1. 



Find Dp so for n > 1, 



+ 1 1 < Dp^n 



2{n + 2^) 

and (TTTTIl implies (Q- □ 

Proof of Theorem \7.1\ By ()A.16|1 . we can suppose d^s = 0. By scaling 
and translation, we can suppose xq = 0, a = 1. By Theorem IA.2| we 
need to get lower bounds on J^-^\iTn{x)\'^w{x) dx. By Holder's inequal- 
ity, for any a, P, p & (1, oo), and q dual to p, 



|7r„,(x)|" dx 



< 



\'7T„(x)\°'w(x)'^w(x) dx 



1 \ 1/p 

Mx)rw{xfp 



wix] 



dx 



We want to pick g, a so /?g = r, ap = 2, Pp = 1, that is, 

1 + r 2r ^ r 



g = 1 + r 



a 



The result is that 
1 

2 



7r„(x)| (ix > I / w{x) dx 



1 



1 + r 

-l/r 



1 + r 



1/9 

(7.12) 
(7.13) 



1 \ (i+r-)/'- 

|vr„(x)|2''/(i+'-) dx 



>C\7inm'n-'~''" (7.14) 
Taking the inf over all 7r„'s with 7r„(0) = 1 and using ()A.5|) . we get 

i^n(0,0)-^ > Cn-^-'''' 
which is (Q. □ 
Example 7.3. Let d/i be the measure on [—1, 1] given by 
dii{x) = Ca,6|x|"(l- IxH^rfx 



where a > 0, 6 > — 1, and C is a normalization constant. 
This is an even measure so P2n-i(0) = 0. Moreover, 

P2n{x) = qn{x^) 



(7.15) 



(7.16) 



32 Y. LAST AND B. SIMON 

where g„ are the OPs for the measure obtained from an x ^ y = x"^ 
change of variables. Since dx = {dy)/y^/'^, we see qn are the orthogonal 
polynomials for the measure 

Ca,b\yr^'''^'i^-yfdy (7.17) 

Thus, up to a constant, 

ir„(0,0;c//i) = ir„,,.^(l,l) 

where Kn;a,i3 is the CD kernel for the Jacobi polynomial associated to 
(1 — x)"{l + xY dx with a = ia — |, /3 = 6. If we call these orthogonal 
polynomials Jq^/j, and Ja^p the conventional normalization, then j5Hl3T] 

\\JaA\2--C^o:ln-^ 



SO 



and 

1) ~ 

Taking a = i(a — 1), we get 

ir„(0,0;d/i) (7.18) 

We can take r in Theorem 17.11 arbitrary with ra < 1, so ()7.H|1 cannot 
be improved. □ 

The following shows that in some cases the power of n in Theorem 17. II 
is optimal: 

Theorem 7.4. Let dfi{x) = w{x) dx where supp((i/i) C [—1, 1] and 

k(a;)|<C|x|" (7.19) 

for some a < 1. Then 

|i^n(0,0)| >Cini+° (7.20) 

Remark. For w{x) = Ca\x\" on [—1, 1], ()7.2|1 holds for any r < 1/a, so 
(l7.2Up says (|7.3p cannot hold for any smaller power of n in case r > 1. 



Proof. Let 7r„ be the polynomial of Theorem IA.7I where xq = 0, a = 1. 
On|x|G[^,^],j = 0,l,2,...,n-l, 



k„(a;)| < 



1 j = 

2n ^ 2j 



' +^ J = l,2,... 
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SO, by (Hini), 



1 1 

+ 



2^2 2f 
<2c|^ + J-fl + Vil±Il: 



n 



since in? > and YlT=i ^^^-^ ^ ^ since a < 1. 
By deg(7r„) = 2n — 2, and Theorem IA.21 

which is (TT^ . □ 

8. Upper Bounds via OP Methods 

Our main purpose in this section is to note that the upper bounds 
produced by the method of Erdos-Turan 19] provide a universal bound. 
So long as d/ig = near xq G supp((i/i) and w{x) is continuous and 
nonvanishing at Xq, the bound is independent of the value of w aX Xq\ 
Upper bounds on spacing imply lower bounds on the density of zeros. 
Deift-Simon 8j obtained universal lower bounds on the density of zeros, 
so the bounds we find are a kind of microscopic analog of theirs. One 
key to the Erdos-Turan method is 

Lemma 8.1. Let Xi < X2 < ■ ■ ■ < Xn in and let 1 < j < n — 1. 

Then there exists a polynomial vr of degree at most n — 1 so 

n{xe) = l<i<n;i=^j,j + l (8.1) 
7r(xj) = 7r(xj+i) = 1 (8.2) 
n{y) > 1 m [xj,Xj+i] (8.3) 



Proof. Let 
If 

take 



M^)= n i^-^j) (8-4) 

TCo{xj) = 7ro(xj+i) (8.5) 



X = —— 8.6 

7ro{xj) 
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SO dHU/ra hold (we will look at dHU shortly). If dHS)) fails, for 

y e M.\[xj,Xj+i\, let 

TCy{x) = {x- y)no{x) (8.7) 

As y runs through {—oo,Xj], runs from 1 down to 0, and as y 

runs from oo to xj+i, the ratio runs from 1 to oo. Since ()8.5j) fails and 
7ro(xj) and no{xj+i) have the same sign, there is a unique y with 

7iy{xj) = ny{xj+i) (8.8) 

so take 

In any event, tt obeys 1)8.111 and ()8.2)1 . By Snell's theorem, vr' has a 
zero between any two zeros of vr, and so by counting degrees, exactly 
1. It follows that tt' has a local maximum in [a;j,a;j+i] and no local 
minimum, so ()8.3p holds. □ 

Let {xj}^^i be the zeros of the OP, P„, associated to a measure dfj,. 
Recall (see Theorem IA.4|) that there are positive weights {Aj}"^^ so 

/n 
n{x)dfi{x) = ''^^ \ef[{xi) (8.10) 

for any polynomial vf with degvf <2n — 1. 

Theorem 8.2. For any j = 1,2, . . . ,n — 1, 

fi{[xj, Xj+i]) < Xj + Xj+i (8.11) 

Proof. Let tt be the polynomial of degree n — 1 or less given by 
Lemma f8. II Let tt = tt^, so degTf <2n — 2<2n — l. Since 7r{xe) = 0, 
^ 7^ J, J + 1, and n{xj) = Tr{xj+i) = 1, 

RHS of (pnnil = Xj + Aj+i (8.12) 

Since vr > and tt > 1 on 

LHS of (IHJTID > fi{[xj, Xj+i\) (8.13) 

so (IHinil implies (IH1T| . □ 

To exploit ()8.11|) . we need upper bounds on Xj. 
Suppose 

i?± = LTsupp(d/x) (8.14) 

and for [E^,E+], 

d{E) = max(E+ - E,E - E^) < E+ - E^ (8.15) 
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Theorem 8.3. Suppose I is a closed interval on which dfi is purely 
a.c. and max^.g/ w{x) = < oo. Then for each 6 > and all weights 
Xj associated with Xj G / and dist(xj, ]R\/) > 6, we have with m = [|] 
and n the number of zeros, 

w^nma^rdjE) ^^^^^^ 
m \n'^ J 

where the C'(^) is uniform in all A 's with the given 5 {and depends on 
max/ d{E)). 

Proof. Let a = ma.xj d{E) and xo be the Xj for Xj, and let 7f(x) = 
Hm{x;xo,a) given by Theorem IA.71 Since 7r(x) = 1 and deg('7f^) < 
2n — 1, we have 

Xj < ^ A^7r(x^)^ 

dfi{x)fc{xy = K1 + K2 

where Ki is the integral over (xq — 6,Xo + 6), and K2 the integral over 
all other x. By (jOH), 

with estimates only dependent on a and 6. For n{x)^ = O(^) on the 
region of integration and //(M) = 1. 

If (xo — 6,xo + 6) C /, d^ < w+ dx, so ()8.1fij) follows from ()A.35j) . □ 

Theorem 8.4. Suppose I is a closed interval on which dfj, is purely 
a.c. and 

< w_ = min w{x) < max w{x) = w+ < 00 (8-17) 
Then for any E G J™*, 
limsup n[zl^\E) - zI~^\E)] < 47rd{E) lim 

miiii «;(X'j I IX — -C/ 1 

(8.18) 

In particular, if E is a point of continuity of w, 

LHS of < A7Td{E) 

independently of the value ofw{E). 
Proof. Clearly, 

\zl^\E) - zt'HE)\ 

< [mm{w{x) I zi-'\E) <x< z^^\E)}]-' f,{zi-'\E) , zI^\E)) 

(8.19) 



max{ti?(x) 




\x-E\< 6} 


min{w(x) 




x-E\ 


<6} 
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From this, (jH.llj) . ()8.1fjj) . and limn/m = 2, we get ()8.18|1 by using 
the fact that since E is a hmit point of an infinity of zeros, we have 
hm„^ookn'(^)-^| =0. □ 

This is the promised universal lower bound on the density of zeros. 
The method is flexible enough to say something if w{x) has a zero of a 
fixed order. 

Theorem 8.5. Suppose dfi is purely absolutely continuous in a neigh- 
borhood of Eq, and for some q > 0, 

< 7_ = liminf < limsup . = 7+ < cx) (8.20) 

' x^Eo \x-Eo\<^- ^^Eo \x-Eo\l ^ ' 

Then 

limsup n|4^)(^o) - 4"'H^o)| < 00 (8.21) 
Proof. By fl8.20j) for any 6, there is so for n > N, 

K[^i-'\Eo),z^:\Eo)]) > {^^-6)2"^{q + ir\zl^HEo)-zi-'\Eo)r' 

(8.22) 

By the proof of ()8.16|) . the A's associated to z^{E) obey 

^ . < Ci(7+ + 6)[maxiC,n-\ |z^(Eo) - Eo\)]'^ 
^ ~ n 
for constants Ci, ^2- (^^ . (^0^ . and imply (jHISH)- □ 

Given our bounds in Section [71 we can also say something when the 
singularity of the weight is not as regular as some power. The key is 
an abstraction of an argument of Nevai ^1] (see also Golinskii 



Theorem 8.6. Let a = max(supp((i/i)) — min(supp(c?yu)). Fix integers 
p, q so that 

{2p - 2)2« < 2n - 1 (8.24) 
Then for any successive zeros E,E' of Pn, we have 

\E -E'\<^ [K,i^ {E + E'), \ (E + E'))]'^'' (8.25) 

Proof. Let vf be defined in terms of the vr of Theorem I A. 71 by Theo- 
rem 

n{x) = [np{x;l{E + E'),a)]'i (8.26) 
By iK7^ . deg[frf < 2n - 1 so, by (jHHUl), 

/n 
mx)\'df^{x) = J2^ME,)\' (8.27) 

.7=1 
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since ^ A^- = 1, > 0, and min\Ej - \{E + E')\ = 1\{E - E')\. 
Since tt{UE + E')) = 1, by Theorem 



KX2iE + E'),'^{E + E'))-'< j\ 



TT 



(8.29) 



Since \E — E'\ < a, 



1 a a 

2^^ 2p\E - E'\ - p\E-E'\ ^' ' 

is immediate from (jHSHI), and ^Ml- □ 

The following abstracts an argument of Golinskii, who needed to 
make global hypotheses since he relied on estimates of Geronimus: 

Corollary 8.7. Suppose that for some interval I, A > and C , we 

have 



suY>\Kn{E,E) \ < C{n + l) 
Eel 



Then for any 6 > 0, 



lim sup \E — E'\ 

E,E' successive zeros 
dist(£;,]I 



n 



logn 



< 00 (8.31) 



Proof. Pick q = [log n] and p as large as possible so that (j8.24j) holds. 
Since 

(n + 1)^ = exp(Alog(n + 1)) 
is bounded and implies (jHZE))- □ 

Combining this corollary and Theorem 17. 11 we obtain a local version 
of Golinskii's [\A^ result: 

Corollary 8.8. // dZH) and (Q hold, then we have (IHlIT|l for I = 

{xo - a,xo + a). 

We also have the following (a local version of results of Nevai 
and Golinskii 



Theorem 8.9. Suppose for some interval I we have that 

dn = wdx + dfis (8.32) 

where 

log wdx>— 00 (8.33) 
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Then for any 5 > 0, 

lim sup \E - E'\n^''^ <oo (8.34) 

"'^'^ E,E' successive zeros 
dist{E,R\I)>& 

We need the following lemma: 
Lemma 8.10. If J has Jacobi parameters obeying 

oo 

^|a„-l|2 + |6„|2<oo (8.35) 

n=l 

then for any S > 0, there is Cs so 



sup \\T4E)\\ < exp{CsV^^) (8.36) 

£■£[-2+5,2-5] 

Proof. Define u^{E) = e*"^ where 2 cos 6' = E and < 6 < n. By 
standard variation of parameters about u^, one proves 

n 

LHS of < Y[{1 + C{\bj\ + \aj - 1|)} 



Since 1 + x < exp(x) and 

n r " 

5](|&,| + |a,-l|)< 5^(|6,f + |a,-l| 

i=i i=i 



1/2 



[2n 



1/2 



(j8.36p is immediate. □ 

Proof of Theorem \8.yi By scaling, we suppose / = [—2,2]. By Corol- 
lary 

Kn{x, x; dfi) < Kn{x, x; du) (8.37) 

where 

du = Xl-2,2]W dx (8.38) 

Let du be the normalized du. By the theorem of Killip-Simon ^25j . 
the Jacobi parameters obey ()8.35|) . so by ()8.37p . 



sup Kn{x, x; dfj,) < exp(^Csy/n + 1 ) 

xe[-2+5,2-5] 

In Theorem 18.61 take q = [v^] P large as can be so ()8.24j) 
holds. Then p ~ c^/n and ()8.25|) implies ()8.34|) . □ 
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9. Lower Bounds via OP Methods 



In this section, we will get lower bounds in terms of the CD kernel 
alone. The basic method is due to Golinskii [T3], but when he applied 
the method, he made global assumptions on the measure, and we want 
to note that local assumptions suffice. Other OP lower bound methods 
are due to Erdos-Turan 9\ and Nevai 



Theorem 9.1. If E,E' are distinct zeros of Pn{x), if E = ^{E + E') 

and 6 > j\E — E'\, then 



\E-E'\ > 



[6^-i^-E'\) 
3n 



Kn{E,E) 



.1/2 



Kniy,y) 



(9.1) 



Remarks. 1. In most applications, 6 is fixed and \E — E'\ 0, so 
(5^ — {^\E — E'l)'^ ~ 5^ > 0. In typical cases, the inf and sup of Kn{y, y) 
for \y — E\ < 5 are comparable and 1)9.11) gives an 1/n lower bound. 

2. This theorem also yields a result with the same asymptotics as 
()2.22|) for Kn{E, E) > 1, while the sup is bounded exponentially in n. 

3. Interestingly enough, the proof here depends on ()2.19|) written as 

4. It is interesting to compare Theorems 19 . II and 12 . 21 ()2.17p only de- 
pends on information at Eq while i\9.1\i has a sup Kn{y, y) over a neigh- 
borhood, but ()2.17|) requires information on both solutions of ()2.H) 
while ()9.ip only on p„. 



Proof Since Pn{E) = Pn{E') = and E ^ E', we have 

Kn{E,E') = 



(9.2) 

by ()A.3|) (the Christoffel-Darboux formula). Thus, (supposing E < E' 
for notational convenience), 



K^{E, E) = K^{E, E) - K^{E, E') 



< \E 

< \E 



E'. 

\y 

E'\{5^ 



sup 

y-E\<hE-E' 



sup [6'^ 

\y-E\<5 



i(E-E')Fr^ 
{y-m'^' 



(9.3) 

(9.4) 

(9.5) 



Kn{E, y) is a polynomial in y of degree n, so by ()A.29|) and ()9.5|) . 

Kn{E,E)<\E-E'\{6^-[\{E-E')f}-\?,n) sup \K^{E,y)\ (9.6) 



\y-E\<S 
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By the Schwartz inequahty, 

sup \K4E,y)\<K,,,{E,EY/' sup \Kr,{y,y)\'/' (9.7) 

\y-E\<5 \y-E\<5 

dnH) plus dnUI) imply dHU). □ 
Corollary 9.2. Let I be an interval on which 

t= sup \\Tr,{y)\\ < oo (9.8) 

n,y€l 



Then, for any E G / , 



r2 



liminf n\z^+^\E) - z\^-^\E)\ > dist(E,M\/)2 (9.9) 

n— >oo 3 

Remark. This should be compared with what follows from ()2.15j) and 
()2.17|1 which implies 

LHS of ^ > t-^ 

Proof. fl9.9p follows from ()9.H) if one notes that for y & I, 

{n + l)t-^ <K„{y,y)<{n + l)t^ □ 

We can also use Theorem 19.11 to get a lower bound in terms of local 
bounds on the weights. 

Theorem 9.3. Suppose dfi = w dx + dfi^ where d^s{xo — 6,Xo + 6) = 
and 

< inf w{x) < sup w{x) < oo (9.10) 

\y-xo\<S |y-a;o|<<5 

Then for any e < 6, 

inf liminf n|4^^)(?/) - 4"^'(y)l > (9.11) 

\y-xo\<e n^oo 

Proof. By fl9.ip . it suffices to prove 

sup [n-^K^{y,y)] <oo (9.12) 

|?;-xo|<e 

and 

inf [n-'K^{y,y)]>0 (9.13) 

\y-xo\<e 

By ()9.10|) . for e fixed, uniformly in y with \y ~ Xq\ < e, we can 
find a fixed scaling and some translate of c(4 — a;^)^/^X[-2.2] dx lying 
below dfi{x — y), so using ()A.14j) and the explicit Kn for Chebyshev 
polynomials of the second kind (i.e., the free K^), we get (j9.12j) . 

On the other hand, by (IA.4I1 / (IA.5I1 . 

Kn{xo,Xo) > (^j |7r„(x)pd/i(x)^ (9.14) 
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for any 7r„ of degree n with 7r„(xo) = 1. Using a suitable 7r„/2 of the 
form given by Theorem lA . 71 and estimates we used earher in this paper, 
we get an 0{n) lower bound on K^, that is, ()9.13|) holds. □ 

10. Zeros of POPUC 

While we have discussed OPRL up to now, virtually all the ideas 
extend to POPUC. POPUC are defined by taking the first n — 1 re- 
cursion parameters (Verblunsky coefficients), o;o,---,«n-2 in D, and 
picking (3 G and letting 

(pn{z) = Z(pn-l{z) - (^^l_i{z) (10.1) 

Since ^Pn_i is nonvanishing on © and Iv^n-il = Iv'n-il on dll>, we have 
\(Pn-i/^Pn-i\ < 1 on D by the maximum principle. Thus, (^„(z) is 
nonvanishing on D and, by symmetry ((^„(l/z) = z~"-[ip'^_-^ — Pzipn-i]), 
nonvanishing on C\D. Thus, the zeros of lie on dH); indeed, they are 
the zeros of a finite unitary matrix (see Theorem 8.2.7 of [H^). Zeros 
of POPUC are discussed extensively in Golinskii [13], Cantero- Moral- 
Velazquez Hj , Simon i46j , and Wong [33] . 

As explained in Section 10.8 of 130] , there is an OPUC analog of 
(I2I7|)-(I211) (namely, (10.8.3)-(10.8.5) of [30]) which immediately leads 
to an analog of Theorem 12. II While (10.8.3)-(10.8.5) are stated for the 
solutions ip. + F(p., they also hold for ip. and if.. Key to this analog 
is the orthogonality of ^p.{^i) ^.{z2) for two zeros of (pn{z)- This 
follows from the CD formula for OPUC (see Theorem 2.2.7 of ^3^) for 
(pn{zi) = (pn{z2) = implies Zjipn{zj) = P(pl{zj), so 

Ki^l) <(^2) - ^12:2 Vn{Zl) <^n{z2) = 

and thus, by (2.2.42) of [39j, 

n 

J2^^^,{Z2)=0 (10.2) 
i=o 

Combined with the techniques of Section 12.1 of [30] and our proof 
of Theorem 13. 4| we get 

Theorem 10.1. Let {a„}^Q be a set of Verblunsky coefficients that 
obeys 

oo 

^ |a„+i - «„! < oo (10.3) 

n=0 

and 

Ittnl^O (10.4) 
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Then the zeros of the POPUC, (pn{z), for any choice of 13 have uniform 
clock behavior on any compact subset of d&\{l}. 

Remarks. 1. An interesting example is an = {n + 2)~^ for any /5 > 0. 
This is related to a conjecture of jUj, albeit the conjecture there is for 
OPUC, not POPUC. 

2. The density of zeros in this case is d6/27i on see Theorem 8.2.7 
and Example 8.2.8 of [39j. 

3. If (fTTD^ / fTinjl are replaced by 

oo 

< oo (10.5) 

n=0 

then c?D\{l} can be replaced by 9D. This is a result of j41j. Because 
we have global control in this case, one does not need a priori 1/n 
bounds on zero spacing. 

There are also analogs of the bounds of Sections El El and El 

1. One has that 

^Me-e«| = -p^ (10,6) 

SO there is a bound like ()5.4|) for POPUC (all of whose zeros lie on 
and thus, there is an analog of Theorem 15.11 

2. If one defines rin{z) by 



then 



This follows from (2.2.71) of which implies 
d 



(10.8) 



^log|(^„+i(re^^)r 



= {n + l) + \^n+,{e^')\-'J2Me^')\' (10-9) 

'•=1 j=0 



By the Cauchy-Riemann equations 

LHS of mM = 2 ^ arg[<^„+i(e*^)] (10.10) 

Since fll0.7j) implies 

r/„ = - (n - 1)0 + 2 arg[(/.„_i(e^^)] (10.11) 

we obtain (fl7r8|l . 
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()1().8|1 implies drjn/dO > and, given that zeros of ipn occur when 
rjn = arg/5 (mod 27r), bounds hke those of Theorem 16 .31 on zero spacing 
for POPUCs. 

3. Since all the techniques of the Appendix extend to OPUC, the 
estimates of Sections |H1 and El extend to POPUC; indeed, somewhat 
weaker variants occur already in Golinskii (14j . 

11. Examples, Counterexamples, Conjectures, and 

Questions 

(a) Clock behavior based only on local behavior. Suppose ()7.1|1 holds 
on for some (c, d) C M, we have dfis{[c, d]) = 0, w > on {c, d) and w 
is C°° there. We have proven 0{l/n) upper and lower bounds in this 
case. This leads to the natural question: 

Open Question 11.1. Under the above hypothesis, does one have 
clock behavior on {c,d)7 

This is a very subtle question because clock behavior involves the 
density of states, and it is not even clear that exists on (c, d) only 
under the above hypothesis. What is clear is that if the density of 
states exists, it is a global quantity and not just dependent on w on 
{c,d). We want to demonstrate this by example. We will need the 
following: 

Proposition 11.2. Let dfiQ be given by 

dfio = (27r)~V4 - x2 X[-2,2] dx 

Suppose f is a C"^ function on [—2, 2] with f'{2) 
n-2) = f"{2), and 

f>a>0 j fdfXo = l 

where a > is a positive real. Let 

dfi{x) = f{x) dfiQ^x) (11-3) 

Then the density of states exists for dfi and is given by p.8|) and there 
is clock behavior uniformly on each interval [—2 + e,2 — e]. 

Proof. By Theorem 13.2.1 of jlU], there is a map Sz2 of real measures, 
dp, on (9D (i.e., those measures with real Verblunsky coefficients) to 
those measures d/i on [—2, 2] which are of the form / dpo where 

j f{x){A-xY^/Ux <oo (11.4) 



(11.1) 

= f (-2) = 0, 

(11.2) 



44 



Y. LAST AND B. SIMON 



and under this map, dp{0) = w{6) dO where 

w{e) = cf {2 co^e) (11.5) 

Thus, w is (including aX 6 = 0, vr) and so Xl^-oo ''^^I'^^nP < 
oo, which imphes ^.^.ool""^"! < Thus, by Baxter's theorem 

(Theorem 5.2.1 of [39j), the Verblunsky coefficients are in . By 
(13.2.20)/(13.2.21) of (irTH|l holds, which implies the claimed re- 

suh. □ 

Example 11.3. Let djiQ be d/ig scaled to [—1, 1]. We pick /i obeying 
the hypothesis of Proposition II 1 . 21 and /2 scaled to [—1,1], and so that 
dill = fi dfio and dfi2 = f2 djlo obey djjii = dfi2 on [— |, |]. Both have 
clock behavior on [— i, i] but with different density of states, namely 
(0|) and (HIHl) scaled. □ 

(b) Pointwise upper bounds. We obtained lower bounds on Zn\Eo) — 
zI^^\Eo) if Tn{Eo) is bounded, but our upper bounds required control 
of Tn{E) for E in a. neighborhood of E^. 

Open Question 11.4. Are there upper bounds on spacing if we only 
know that T„(i?o) is bounded? 

(c) Improved spacing estimates. 

Open Question 11.5. Can ^/n in ()2.25|) be improved? 

(d) More on spacing and '~i{E). We saw that ^{Eq) = plus some 
regularity of u at Eq implies an 0{l/n) lower bound. Does it imply 
clock spacing? In particular. 

Open Question 11.6. Is there local clock behavior for a.e. Eq with 
7(£'o) = in the case of almost periodic Jacobi parameters? 

Example 11.7. |2Z| has proven, for a < |, the OPUC analog of 
Poisson behavior for a„ = 1, 6„ independent random variable of the 
form bn = Cn~'^Wn where Wn is uniformly distributed in [—1,1]. We 
assume their result is true in the Jacobi case. Of course, j{E) = in 
this case. We do not have a contradiction with Corollarv 14. 21 since this 
model is not ergodic. The example does show though that ergodicity 
is a necessary hypothesis. □ 

Example 11.8. Corollarv l4 . 21 shows that ergodicity along with Poisson 
behavior imply positive Lyapunov exponent. This raises the natural 
question: Does ergodicity along with a positive Lyapunov exponent 
imply Poisson behavior? The answer is negative, as can be shown by 
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the following example: Consider the Jacobi matrix with = 1, 6„ = 
A cos(27ra;n+^), where |A| > 2 and a is a (Liouville) irrational for which 
there is a sequence of rationals {pn/ <in}'^=i such that |«— Pn/*?™! < 
This is an ergodic Jacobi matrix and it is well known (see, e.g., [IJ that 
its Lyapunov exponent 7(-E) is positive for any E. By using the results 
of Avron-van Mouche-Simon [2, and considering scales of the form mg„, 
where m > 2 is an integer, one can show that, for each 9, there would 
be clusters of m — 1 zeros each of which is contained in an interval 
whose length is of order (2/|A|)'^"''^. As is varied, these clusters will 
move over regions whose size is roughly of order This behavior 

contradicts Poisson behavior. More precisely, it is possible to show that 
Poisson behavior does not occur for Lebesgue a.e. E in the spectrum. 

(e) Zero spacing and the Szego condition. 

Open Question 11.9. Does one have 0{n~^) bounds (upper and 
lower) when a Szego or quasi-Szego condition holds? 

(f) Spacing at zeros of w{x). 

Open Question 11.10. What can one say at zero spacing at points 
xo where w{x) has a "regular" zero, that is, w{x) ~ |a; — xo|" for some 
a > 0? 

(g) Edge zeros when a„ = 1 — rT'^ . The following illuminates Theo- 
rem EUl 

Example 11.11. Let = 0, = 1 — n'^"' for 7 > 0. Then Theo- 
rem EiU applies and there is clock behavior away from —2 and 2. If 7 > 
2, Theorem 13.81 applies and the largest Ej has Ej = 2 — Cn""^ + o(n~^). 
We claim for 7 in general 

2 - C2n-"' + 0{n-^) < Ef"" < 2 - Cm"^ (11.6) 

capturing the leading behavior for 7 G (0,1). The upper bound in 
flll.6|) comes from monotonicity of E^^^ in the a's and the fact that 
for J„;i? maxi<j<„ I Oj I ~ 1 — The lower bound comes with a trial 
vector that lives in [n/2, 2] with maximum at 3n/4 and constant slope 
in between. □ 

(h) OPUC. This paper has a fairly complete analysis of OPRL and 
POPUC. Many questions remain for general OPUC. 
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Appendix A. Tools of the OP Trade 

As explained in the introduction, this paper is intended for two au- 
diences, so we include here a summary of tools well known to the 
OP community but not so well to the Schrodinger operator commu- 
nity. Because the tools, while powerful, are simple, we can even give 
complete proofs. We will discuss the Christoffel variational principle. 
Gauss- Jacobi quadrature, Bernstein's inequality, and Dirichlet-Fejer 
trial polynomials. 

(a) Christoffel variational principle. We will define OPRL for 
arbitrary positive measures (with finite moments) even if /^(IR) 7^ 1. 
The monic polynomials, Pn, are independent of normalization, but 
the orthonormal polynomials, p„, are not. For example, po{x, dfj,) = 
dfi(R)-^/^. The Christoffel- Darhoux (a.k.a. CD) kernel or reproducing 
kernel is defined by 

n 

Kn{x,y) = pj{y)pj{x) (A.l) 
i=o 

We will use Kn{x, y; dfi) if the measure needs to be more explicit. The 
name reproducing kernel comes from 

(P„/)(x) = j K„{x,y)fiy)dM (A.2) 

where P„ is the projection in L^(M, dfi) onto the space of polynomials 
of degree n. 

We need the following in Section IHl 



Theorem A.l (CD formula). We have 
Kn{x,y) = a„+i 



Pn+l{y)Pn{x) - Pn{y) Pn+l{x) 



y~x 



(A.3) 



Proof. This is a discrete version of integrating a Wronskian. Take the 
equation ()1.3p at y, multiply by Pnix) and subtract ()1.3|) at x, multi- 
plied by Pn{y), and obtain 



Qn+lix,y) =Pniy)Pnix) + Qn{.X,y) 

where Qn+i is the right side of ()A.3|) . Since p-i(x) = 0, Qo = 0, so 
()A.3|) follows by iteration. □ 

Remark. For x, y real, pn{x) , pn{y) , are real so the bars are not needed. 
Indeed, one can drop all the bars for complex x, y, but given the OPUC 
analogs, it is natural to use the bars. 
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Xn{xo;dfi) = inf( / |7rn(a;)|^ 
Then 



deg7r„ < ra; 7r„(xo) = 1 (A.4) 



A„(xo; djj,) = Kn{xo, Xq; dfi) ^ 



(A.5) 



Remarks. 1. A„ are called Christoffel numbers. More generally, we have 
p-Christoffel numbers defined, for < p < oo, by 



A„(xo,p;rf/i) = inf ( \nn{x)\P dfi{x) 



deg7r„ < n; Hn{xo) = 1 (A.6) 



2. Another way of writing ()A.4|) is that A„ is the optimal constant 



m 



or 



\TTn{Xo)\^ < Xn{Xo;dfl) ^ / \'Kn{x)\^ dfl{x 



|7r„(xo)|^ < A„,(xo,p; rfyu) ^ / |7r„(x) |p (i/i(a; 



(A.7) 



3. Our proof shows the inf in ()A.4|1 is a min and the minimizing vr 
is given by 

7l{x) = K{x, Xq) 



Proof. Expand 7r„ in terms of the orthonormal basis {pj}"^o- 



x] 



j=0 

T!'n{xQ) = 1 is equivalent to 

n 

^ajPjixo) = 1 
By the Schwartz inequality 



(A9) 



(A.IO) 



1 < K{xo,xo) ^ 

j=0 



a] 



K{Xo,Xo) / 7ln{xYdfl{x) 



(All) 



(A.12) 



(A.13) 



where equality occurs in ()A.12|) if aj = pj{xo)/K{xo,Xo), that is, if 7r„ 
is given by ()A.10|) . ()A.5|) is immediate from this case of equality and 
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Christoffel numbers have been a critical tool in OP theory for over 
a century, with important uses by Erdos-Turan [S] and turned to high 
art by Freud and Nevai; see Nevai They can be used for lower 

bounds on K, that is, upper bounds for A„, by using any convenient 
trial polynomial for 7r„ (see (d) below). One gets upper bounds on K, 
that is, lower bounds for A„, by the immediate 

Corollary A. 3. // dfx > du, then 

Kn{xo, Xq; dfi) < Kn{xo, Xq; du) (A. 14) 

Remarks. 1. This shows the true power of Theorem IA.2I and the need 
to allow /i(M) 7^ 1. 
2. In particular, if 

dfj, = f{x) dx + djjis (A. 15) 

then 

i^„(xo, xo; dii) < K„{xo, Xq; f dx) (A. 16) 

(b) Gauss-Jacobi quadrature. The main result here is 

Theorem A. 4 (Gauss-Jacobi quadrature). Let /i be an arbitrary pos- 
itive nontrivial measure on M with finite moments. Fix n and define 
dfin to be the point measure with weights only at the zeros {x^"''}"^]^ of 
Pn{x) and weights 

dfini{x^;^}) = Xnix^;^;dfi) (A.17) 

the Christoffel numbers of dfi. Then, if tt is a polynomial of degree 
2n — 1 or less, we have 

J 7c{x) dfi{x) = J n(x) dfin{x) (A.18) 

Remark. In our applications, we will care much more that the masses 
are at the zeros than the variational formulae for the weights. 

Sketch of Proof. Here is a proof intended for Schrodinger operator ex- 
perts. (For the more usual OP proof, see Freud's book [IHl-) Let Jn-F 
be an X matrix in the upper corner of the Jacobi matrix, (jl.4p . 
associated to dfi. Then the recursion ()1.3p implies that if 

Uj{z) = p,^i{z) j = l,2,...,n (A.19) 

then 

[{Jn;F - Z)u]j = -an+lSjnPn{z) (A. 20) 

Thus the eigenvalues of Jn-p are the zeros of p„ and the normalized 
eigenvectors are pj-i{z)/ K{z, zY^"^ . It follows that djjn is the spectral 
measure of Jn-p with eigenvector (1 ... 0)*. 
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For any measure drj and associated Jacobi matrix, J, since dr] is the 
spectral measure for 5 = (1 ...)*, we have 

xUr]={5,f5) (A.21) 

so / x^^^ drj = {J^5,J^5) for = 0, ...,n — 1 depends only on 
{J^5}1'1q and so on {afc,6fc}^~Q. Moreover, 

d7] = (A.22) 

only depends on {ofc, &fc}fc=o U {^n}- Thus, J and J„;f, which have the 
same set of these parameters, have the same moments of order up to 
2n — 1, that is, 

j x^d^i = j x^ dfXn-^F 0<k<2n-l (A.23) 

which is ^EdB- n 

(c) Bernstein's inequality. These inequalities control tt^ in terms 
of n and 7r„ for polynomials 7r„ of degree at most n. 

Theorem A. 5. Let d3 be the unit circle in C, d3 = {z \ \z\ = 1}. 

Let Tin he a polynomial of degree n. Then 

sup |7r^(-2)| < n sup |vr„(2;)| (A. 24) 

Remark. If TTn{z) = z^, one has equality in ()A.24|) . 
Proof. (Szego [20]) Since vr„(e*^) = YTj=o'^j^^^^ ^ have 



so 



I-21T " 



ZTT 



F^{e-ip)e'-^'-^>n^^ie'n^ (A.26) 

In 



where 



n-l 

F„(0)= 5^ {n-\j\)e^^' (A.27) 

j = -71+l 

(for the j > terms in ()A.26j) . integrate to zero). 
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By cancellation, 

(1 - cos e)Fn{e) = 1 - cos{n9) (A.28) 
so Fn{e) > and, by ^K27\ . J = n. Thus, 

K(e*')| < Iknllooy = n||7r„||oo □ 

Theorem A. 6. Let 7r„ be an arbitrary polynomial of degree n. Then 
sup [\'K'^{x)\{a^ -x^f'^]<?,n sup |7r„(a;)| (A.29) 

a,a] x(i[~a,a\ 

Proof. By scaling, we need only check the case a = 2. Define 

%^[Z)= z''T:n{z+-\ (A.30) 



z ^ 

Tin is a polynomial of degree 2n so, by ()A.24|) 



sup K(e'^)| < 2n sup |7r„(a;)| (A.31) 



since e'^ e'^ + e"*^ = 2 cos 6* maps (9D to [-2, 2]. 
By (ESOl), 



SO, by (IXaTll . 



sup |7r^(2 cos6')2 sin6'| < 3n||7r„| 



which is (f09|l for a = 2. □ 

(d) Dirichlet trial polynomials. For use in both (jA31) and (fAT8|l . 

we want a rich set of trial polynomials, 7r„(x). In particular, we want 
7r„'s concentrated near x = xq and otherwise small in some interval 
[xq — a, + a]. By scaling, we may as well consider xq = 0, a = 1. An 
analyst might try (1 — x^)", but that has width n"^/^ — and we will 
see that one can do better. We will get width . One can't do better 
than this, by Bernstein's inequality, if 7r„,(^) = 1 and ||7r„||oo = 1, then 
7r„(a;)>ifor |x|<^-0(^). 

Our choice is related to Dirichlet and Fejer kernels and is, in fact, 
essentially the minimizer for the Christoffel problem with xq = and 
(i/i = ~" x^)'^^"^ dx. 

Theorem A. 7. For any xq G M and a > there exist, for each n, 
polynomials iTnix; Xq, a) so that 
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(i) 
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deg 7r„ = 2n - 2 (A.32) 
vr„(a;o) = 1 (A.33) 



|7r„(x)| < min( 1, 77- + o i " 1) if \x - Xo\ < a (A. 34) 

\ 2n 2n\x — Xo\ J 

(iv) For any 6 < a, 

r^'Mx)\^dx = — + o(^) (A.35) 

Remark. flA.34|) implies |vr„(x)| < Cs/n if |x — xo| > S, and for any e, 
|7r„(a;)| < e if \x — xo\ > Cs/n. 

Proof. By scaling, we can suppose that Xq = 0, a = 1, in which case 
we will call the polynomials D„, that is, 

— a — ] (A. 36) 

Recall there are polynomials T„(a;) (Chebyshev of the first kind) with 

deg T„ = n (A.37) 

so that 

T„(cos^) = cos(n^) (A.38) 

Define Dn by 

^ n— 1 

Dn{x) = - V(-l)%(x) (A.39) 

i=o 

By |X'37j) . (f02|l holds for D^. By fr08|) (and cos0 = ^ ^ = f 
mod tt), T2j(0) = (—1)"', so Dn obeys ()A.33|) for xq = 0. 
By flX3Hll . |T„(x)| < 1 on [-1, 1], so 

\Dnix)\ < 1 on [-1,1] 

which is half of ()A.34|) . For the other half, sum the geometric series to 
see that 

1 (-1)""^ cos((2n - 1)^) , , , 

which implies the other half of ()A.34|1 . 
Since 

r de ^fk = j = o 

/ cos(2A:0) cos(2j^) — = <| ifA;=j^O 

ifMj 
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we have 
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d9 

2^ 



l + -(n-l) 



2n ~'~ 

Since dO = {1 + 0{x^)) dx near x = and 
X = 0, we obtain ()A.35|) when a = 1. 



O(^) away from 
□ 
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